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ABSTRACT 

We simulate the buildup of galaxies by spherical gas accretion through dark-matter 
haloes, subject to the development of virial shocks. We find that a uniform cosmological 
accretion rate turns into a rapidly varying disc buildup rate. The generic sequence 
of events (Shocked- Accretion Massive Burst & Shutdown: SAMBA) consists of four 
distinct phases: (a) continuous cold accretion while the halo is below a threshold mass 
M s h ~ 1O 12 M , (b) tentative quenching of gas supply for ~ 2 Gyr, starting abruptly 
once the halo is ~ M s h and growing a rapidly expanding shock, (c) a massive burst 
due to the big crunch of ~ 10 11 M Q gas in ~ 0.5 Gyr, when the accumulated heated 
gas cools and joins new infalling gas, and (d) a long-term shutdown, enhanced by a 
temporary shock instability in late SAMBAs, those that quench at z ~ 2, burst at 
z~ 1 and end up quenched in 10 12 ~ 13 Af Q haloes today. The quenching and bursting 
occur at all redshifts in galaxies of baryonic mass ~ 10 M@ and involve a substantial 
fraction of this mass. They arise from rather smooth accretion, or minor mergers, 
which, unlike major mergers, may leave the disc intact while being built in a rapid 
pace. The early bursts match observed maximum starbursting discs at z £ 2, predicted 
to reside in 10 13 M Q haloes. The late bursts resemble discy LIRGs at z 1. On the 
other hand, the tentative quenching gives rise to a substantial population of ~ lO n M0 
galaxies with a strongly suppressed star- formation rate at z~ 1-3. The predicted long- 
term shutdown leads to red & dead galaxies in groups. A complete shutdown in more 
massive clusters requires an additional quenching mechanism, as may be provided by 
clumpy accretion. Alternatively, the SAMBA bursts may trigger the AGN activity 
that couples to the hot gas above M s h and helps the required quenching. The SAMBA 
phenomenon is yet to be investigated using cosmological simulations. 

Key words: shock waves — accretion — galaxies: evolution — galaxies: formation 
— galaxies: haloes — dark matter 



1 INTRODUCTION 

Observations reveal a puzzling phenomenon of "max- 
imum starbursting" in massive galaxies at high red- 
shift in cases where major merg ers are ruled out. Ex- 
amples include LIRGs at z ,$ 1 (jHammer et al.l I2005D 
and maximum bursters at z XL 2 (|Genzel et al.l 120061 : 
iForster Schreiber et al.1 [20061 ). In these cases, galaxies 
of stellar masses ~ 10 n M© seem to be forming a large 
fraction of their stars at a star-formation rate SFR~ 100- 
200 MqVT -1 within ~ 0.5 Gyr. This is much shorter than 
the age of the universe then, indicating coherent star 
formation from ~10 11 M Q of gas in the latest stages of 
its collapse/ assembly rather than in its earlier smaller 
progenitors. Furthermore, the burst duration is signifi- 



cantly shorter than the time implied by the typical cos- 
mological accretion rate onto haloes of ~ 10 12 M Q at the 
relevant epochs, and is as short as the time for streaming 
at the virial velocity from the virial radius to the cen- 
ter. Gaseous major mergers, which could have provided 
such rapid bursts, are unlikely in many of these cases, 
where the galaxies are detected to be relatively smooth, 
rotating, thick discs that could not have survived a ma- 
jor merger and show no trace of such an event. At the 
low end of the SFR distribution, most galaxies of stel- 
lar mass Md ^ 10 n M Q tend to be quenched red & 
dead galaxies, and a substantial fraction of such galax- 
ies show strongly suppressed star-formation rates (SFR) 
already at high redshifts, 0.5 < z < 2.7 (jKriek et al.l 
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2006; iNoeske et aTll2007l) . We seek an explanation for 
these phenomena at the two extreme regimes of SFR. 

We report here on new insight gained from spherical 
modeling of gas accretion onto massive galaxies, which 
could be the basis for understanding the above phenom- 
ena and other central issues in galaxy formation. 

Following the classic work on the interplay be- 
tween co oling and dynamical ti mes in ga la xy for- 



mation feees fe Ostrikeil Il977t ISilkl 119771; iBinnevi 
fl977t iWhite fc Reeslll97a iBlumenthal et al.lll984D TwB 

used analytic calculations and hydrodynamical sim- 
ulations to study the evolut ion of a virial shock 
in a spherical configura tion (|Birnboim fc Dekel 120031 : 
iDekel fe Birnboimll2006l ). We predicted the existence of 
a threshold halo mass for the presence of a virial shock, 
at M s h~10 12 M Q , roughly independent of redshift. Less 
massive haloes do not permit a stable shock, as rapid ra- 
diative cooling prevents the hypothetical post-shock gas 
pressure from supporting the shock against global grav- 
itational collapse into the halo center. The gas accret- 
ing through the virial radius of such haloes flows cold 
(~ 10 4 K) into the inner halo, where it may eventually 
shock, feed a disc and efficiently form stars. Once the 
halo grows above M s h, a stable shock emerges from the 
inner halo and rapidly propagates toward the virial ra- 
dius, halting the infalling gas and creating a hot medium 
in quasi-static equilibrium at the halo virial tempera- 
ture. We found that the transition to stability occurs 
when the standard radiative cooling time at the as- 
sumed metallicity equals the time for compression be- 
hind the shock, t comp = (21/5)p/p ~ (4/3)R/V, where 
p is the gas density behind the shock, R is the shock 
radius and V is the infall velocity into the shock. 

Based on these findings, it has been proposed that 
the presence of hot, dilute gas allows the suppres- 
sion of gas supply to the disc, possibly assisting AGN 
feed back, and thus leads to quenching of star forma- 
tion (|Binnevil200l IDekel fc Birnboimll2006f) . This mass 
threshold has proven to be the key for understanding 
the robust division of galaxies into blue, star-forming 
discs v ersus red fc dead spheroid s and their basic prop- 
erties (IDekel fc Birnboirnl 120061 : ICattaneo et~aTl 120061: 
ICroton et al.ll2006h . 

The success of the spherical model in predicting the 
threshold m ass as seen in cosmological hydrodynamica l 
simulations (|Keres et al.ll2005HDekel fc Birnboimll2006l) 
is demonstrated in FigJTJ which summarizes the results 
of iBirnboim et all (|2007t ). It shows as a function of halo 
mass the fraction of cold gas within the virial radius and 
outside the disc in dark-matter haloes from c osmologi- 
cal simulations by A. Kravtsov (described in IKravtsov] 
120031: IKravtsov fc Gned"in] [2005. where they were used 
for other purposes). At all redshifts, there is a transition 
from cold-dominated haloes to hot-dominated haloes 
near 5 x 10 11 M Q . This is compatible with th e threshold 
mass predicted by IDekel fc Birnboirnl ((2006) for realis- 
tic halo metallicities that gradually grow to Z ~ 0.1 
solar today. 

This success of the spherical model motivates our 
present, more detailed and daring investigation of the 
spherical configuration. We study how the development 
of the virial shock converts an otherwise uniform ac- 
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Figure 1. The threshold halo mass for virial shock heating in 
cosmological hydro si mulations, at different redshifts, based on 
(Birnboim et al. 2007). Each point represents the fraction of cold 
gas within 0.1R-R in a different halo, for haloes above a mini- 
mum mass due to numerical resolution. The smaller haloes (open 
symbols) are drawn from a high-resolution simulation in a box 
of comoving size 6/i _1 Mpc, which was terminated at z = 4 and 
therefore do not show counterparts at later redshifts. The more 
massive haloes (solid symbols) come from a simulation of a lower 
resolution in a box of 80 h~ 1 Mpc. Cold gas is defined here by 
T < 10 4 + 0.4(T vir -10 4 ), e.g., T < 2.5 x 10 5 K for M ~ 10 12 M Q , 
where T v ; r ~ 6 X 10 6 K. The transition from cold dominance to 
hot dominance occurs near 5 X 1O 11 M0 at all redshifts (vertical 
line), as predicted by the spherical model for similar metallicities 
jPekel fc Birnboimll2006]) . 



cretion through the virial radius into a rapidly vary- 
ing accretion rate onto the inner disc. The resultant 
disc buildup, and the induced SFR, evolve through 
a generic sequence of events involving two quenching 
episodes and a massive burst, which always starts at the 
crossing of the halo threshold mass, namely in galax- 
ies of Md ~ 10 n MQ at any redshift. We term this 
phenomenon SAMBA, for Shocked- Accretion Massive 
Burst and Shutdown. 

In <J5]we describe the simulation method. In §3\ we 
present the four phases of the SAMBA phenomenon. In 
$4]we dare a tentative comparison to observations. In <J5] 
we describe SAMBA recipes for semi-analytic modeling. 
In SjHl we discuss our results. In related to we 
provide an EPS estimate of the average virial accretion 
rate onto dark haloes. In $B] related to we specify 
the SAMBA recipes for SAMs. 



2 METHOD OF SIMULATIONS 

2.1 The Spherical hydrodynamical code 

Our accurate 1-D Lagrangian hy drodynamical code 
(base d on the code described in IBirnboim fc Dekel 
2003) simulates a spherical gravitating system con- 
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sisting of dark matter and a fraction fb of gas with 
a constant metallicity Z. The initial smooth density- 
perturbation profile is constructed such that it pro- 
duces the desired average accretion rate at the virial 
radius, ^2.21 The dissipationless dark-matter shells de- 
tach from the cosmological expansion, collapse and os- 
cillate into virial equilibrium such that they deepen the 
potential well attracting the dissipating gas shells. The 
gas is cooling radiati vely based on the atomic co oling 
function computed bv lSutherland fc Dopital {1993) , and 
contracts dissipatively into the inner halo. The collapse 
of each gas shell is stopped at ~ 5% of the virial ra- 
dius by an artificial centrifugal force which mimics the 
formation of a central disc. 

The main improvements to th e original code de- 
scribed in lBirnboim fc Dekell ((2003) are as follows. The 
numerical scheme for time evolution now computes the 
pressure and energy of the gas shells within the forth- 
order Runge Kutta solution rather than externally. A 
force corresponding to the cosmological constant A has 
been added to the force equations, making the cal- 
culation fully consistent with the ACDM cosmology, 
though with only minor effects on the evolution inside 
the haloes. The simulations described here consist of 
1,000 gas shells and 5,000 dark-matter shells, corre- 
sponding to a spatial resolution of ~ 50 pc in the inner 
halo growing to ~ 1 kpc near the virial radius. The en- 
ergy conservation level, set by the timesteps and spa- 
tial resolution, is better than 2% over today's Hubble's 
time. Convergence tests of some of the runs with twice 
the spatial resolution yielded very similar results. The 
baryonic fraction /b is assumed to be 10% throughout 
the simulations presented here, crudely taking into ac- 
count certain gas mass loss due to feedback effects. Sim- 
ulations with fb = 0.05% and 0.15% yield consistent 
results. 



2.2 Cosmological accretion rate 

The average fractional accretion rate onto haloes of 
mass M at time t is estimated in §K\ us ing the EPS for- 
malism, following Neis tein et al.l (|2006f ). Eq. (|A5[) reads 



— (M,t) = s(M) — — . 
M y 1 ' y 'DD 



s(M) = 



2/tt 



a 2 (M/q)-a 2 (M) 



1/2 



(1) 

(2) 



Here a(M) is the rms linear density fluctuation on the 
scale corresponding to mass M, linearly extrapolated to 
z = using the linear growth function D(t) normalized 
to unity at to- The constants are 8 C ~ 1.68 from the 
spherical collapse model, and q ~ 2.2 reflecting an in- 
trinsic uncertainty in the EPS formalism. The accretion 
time (M/M) -1 is shown later on in Fig. [5] for haloes of 
10 12 M Q in the ACDM cosmology. A practical approxi- 
mation for haloes of ~ 10 12 M m in ACDM is 



M ( M 

-(M,t) = 0.04 Gyr 



0.15 



The mass growth history of the main progenitor of 
a halo that ends up with mass Mq at time to is given 
by eq. |A"9|) : 

(4) 



M(t\M ) = F^MD- 1 - 1) +F(M )} , 
where J 71-1 is the inverse function of 
dm 



F(M) = f 

J A 



(5) 



M ms(m) 

We specify in $5] the details of this calculation in 
ACDM. 

The initial conditions of the simulation specify 
the mean density profile of the spherical perturbation 
Si(M), referring to spheres of mass M at an initial time 
ti in the linear regime. In order to end up with a halo 
mass Mq at to, and have the proper average virial accre- 
tion rate throughout the halo history, we invert M(t) 
from eq. ^ to a collapse time for mass M, t c (M), and 
then derive the initial profile from 



6i(M) = S c 



D{U 



D[t c {M)] 



(6) 



(1 



(3) 



3 FOUR PHASES OF ACCRETION 
3.1 Generic SAMBA 

Figures [2l [3] and 0] present the SAMBA phenomenon in 
three simulations, spanning a range of final halo masses. 
They display the radii of gas shells as they fall through 
the virial radius and eventually accrete onto the inner 
disc. One may tentatively and very crudely identify the 
disc buildup rate with maximum SFR. We see that ev- 
ery SAMBA consists of four distinct phases, as follows: 

1. Cold accretion. As long as the halo mass is below 
M s h there is no stable shock in the halo. The gas flows 
cold and unperturbed toward the disc following the cos- 
mological accretion rate at the virial radius, which is 
constructed here to be rather uniform. The resultant 
disc accretion rate is growing gradually with time, while 
the specific M^/M^ is declining. For example, in the late 
SAMBA shown in Fig. [21 the cold accretion phase lasts 
till t ~ 4.6 Gyr, with the temperature map showing no 
heating and the disc accretion rate rather uniform. 

2. Tentative quenching. A shock forms near the disc 
upon crossing Af s h, and it propagates outward very 
rapidly. The gas that falls through the shock during its 
rapid expansion is heated and halted to a stop (some- 
times even pulled back) before it cools and rains down. 
A first phase of suppressed disc accretion lasts for 1- 
2 Gyr, e.g. between t ~ 2.7 and 4.3 Gyr in Fig. H 

3. Massive burst. The gas that was slowed down in 
the previous quenching phase now joins new accreting 
gas in a big crunch onto the disc. The accumulation 
of shells in a massive crunch results from the longer 
delay suffered by gas that falls in through the shock 
earlier, when the shock expands faster. The crunch is 
enhanced further when the shock develops a temporary 
instability causing it to contract back to sub- virial radii. 
Thus, ~ 10 11 M Q of gas enter the disc during ~ 0.5- 
1 Gyr, producing a burst of M d ~ 100 - 200Af©yr _1 , 
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Figure 2. Spherical accretion onto the outer dark halo and 
the inner "disc" . Top: Infalling gas shells, with the temperature 
highlighting the evolution of the virial shock. Bottom: baryonic 
mass growth, accretion rate and specific accretion rate (smoothed 
over a virial crossing time), at the virial radius (red curve) and 
at the "disc" (blue curve). A late SAMBA: today's halo mass 
Mo = 4 X 10 12 A/q, metallicity Z = 0.3 and a reduced accre- 
tion rate, leading to a tentative quenching onset at z\ ~ 1.1, a 
burst at zj, ~ 0.7 with efficiency B ~ 5, followed by very effective 
long-term quenching. 

or Md/Md ~ 1 — 2, typically twice the specific virial 
accretion rate. 

4. Post-burst slowdown or shutdown. The halo 
mass is safely above the threshold, with the cooling rate 
significantly slower than the compression rate, allowing 
the shock to join the gradual expansion of the virial 
radius. Bursts that occur prior to z ~ 1.4 are followed 
by a long period of gradually declining specific accretion 
rate, which is comparable to the specific virial accretion 
rate. Bursts that occur after 1.4 are followed by an 
effective long-term shutdown from z~0.7 and on, last- 
ing for several Gyr. The infall is halted effectively while 
crossing the shock during its second expansion phase. 

The timing of the SAMBA events scales with the 
onset time, t\. A relevant characteristic time scale at 
that epoch is the virial crossing time, the time it takes to 
stream at the virial velocity from the halo virial radius 
to the center. Based on the spherical collapse model, it 
is independent of mass and given bjQ 

1 For a proper virial radius, (1 + z) -1 should be replaced by 
(A 2 oo^m0.3^o 7)~ 1 ^ 3 ( 1 + z )^ 1 ■ In a flat universe A (a) ~ [187r 2 - 
82C A (a) - 39f2 A (a) 2 ]/r2 m (a), with a = (1 + z) _1 , Q m (a) given 
in CAl and f2yv = 1 — Q m - 



Red shift: 




Time [Gyr] 

Figure 3. Same as Fig. [2]for an intermediate SAMBA: Mq = 
10 13 Mq, Z = 0.3, typical accretion rate, z\ ~ 2.5, Zb — 1-4 with 
B ~ 6, followed by effective long-term quenching. 



Redshif t 




Time [Gyr] 

Figure 4. Same as Fig. [2] for an early SAMBA: Mq = 10 14 M Q , 
Z = 0.1, typical accretion rate, z\ ~ 4.8, Zb — 2.5 with B ~ 2.3, 
followed by ineffective quenching. 
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^~2.49Gyr(l + z)- 3 / 2 

Vv 



(7) 



In the EdS phase of cosmological evolution (valid to a 
good approximation prior to z~l), i V ir — 0.18tHubbie- 

We define the efficiency of disc accretion by the ra- 
tio of specific accretion rate into the disc and into the 
virial radius, 



B 



Ma/Ma 
M/M 



(8) 



A smoothed version of this quantity over a virial cross- 
ing time, B Vl which can be read as the ratio of the two 
curves shown at the bottom panels of Figs.[2]|U serves us 
for an automatic identification of the four phases. This 
quantity is about unity during the cold-flow phase. We 
mark by t\ the time when B v starts its first drop below 
unity due to the shock onset. Then ti and t% mark the 
subsequent crossings of unity upward and downward, 
which we define as the beginning and the end of the 
burst. Finally, we identify the following upward cross- 
ing of unity, at t±, as the end of the long-term quenching 
phase, when such a phase exists. 

There are several robust SAMBA characteristics 
that only weakly depend on the redshift of the event. 
The onset of the SAMBA is always at the crossing of 
the threshold halo mass, M s h^lO 12 M , thus involving 
a stellar mass Md ~ 10 Mq. The burst peaks about 
half a Hubble time (1.5 — 2.5 Gyr) later, and involves 
,$ 10 11 M Q of rapidly accreting gas in about a virial 
crossing time i V ir (~ 0.5-1 Gyr). The average burst effi- 
ciency is B v ~ 2, with a peak value of _B~6. 

3.2 Systematic redshift dependence 

By simulating the accretion histories of haloes that end 
up with different masses Mq today, we obtain a sequence 
of SAMBAs with a range of starting times, t\, given 
that the SAMBA onset is determined by the crossing of 
a fixed halo mass threshold. Our fiducial case assumes 
the average virial accretion history onto haloes Mq and 
a constant metallicity Z = 0.3 solar; we obtain some- 
what earlier (later) onset times by lowering (raising) 
Z or by raising (lowering) the virial accretion rate at 
the relevant epoch. It turns out that our current fidu- 
cial choice of Z = 0.3 at all times does not follow the 
gradual cos mological growth of metalli city with time as 
estimated in iDekel fc Birnboiml (|2006( ) . This results in 
an artificial overestimate of M s h at high z, by a factor 
~ 2, which we should scale out when considering the 
cosmological redshift dependence of the SAMBAs. 

The characteristic times of 15 simulated SAMBAs 
are shown in Fig. [5] These times naturally scale with 
the Hubble time. The burst duration is always At\, ~ 
0.18t2 — t v i T . The duration of the tentative quench- 
ing is At q ~ 0.8£i for z\ > 2.5 and At q ~ 0.36ti 
for z\ ^ 2.5. The systematic variations with the onset 
time are primarily driven by the gradual decline with 
time of the specific virial accretion rate, which induces 
a drastic change between SAMBAs that start prior to 
z ~ 2.5 and those that start later. This is a distinc- 
tion between the central galaxies of today's clusters of 



B 




Figure 5. Characteristic times of simulated SAMBAs as a func- 
tion of the onset redshift z±, when M = M s ^. The Hubble time 
at the onset of tentative quenching is t\ (open red triangles). The 
ratio tilt\ (filled blue squares) refers to the duration of the tenta- 
tive quenching. The ratio tzjti (filled green circles) refers to the 
duration of the burst, with to, the beginning of long-term quench- 
ing in late SAMBAs. The vertical line marks the transition from 
early to late SAMBAs. The horizontal lines mark the fits used for 
the SAM recipe in §E\ 



> 1O 13 M haloes and the more isolated galaxies in to- 
day's groups of ^ 10 Mq haloes. The former burst 
before z~1.4, when M/M > lGyr -1 , while the latter 
burst later, when M/M < lGyr -1 . In the late SAM- 
BAs, the rapidly expanding shock overshoots to out- 
side the virial radius, where the low accretion rate fails 
to provide enough post-shock compression for support 
against contraction. The resultant temporary instability 
makes the shock drop to sub- virial radii before the pres- 
sure behind it can push it back to the virial radius. Ap- 
parently, this instability appears once the specific virial 
accretion rate drops below 1 Gyr - . 

The effect of the virial accretion rate on the tenta- 
tive shock instability is clearly demonstrated in Figs. [2]- 
0J The strength of the instability as measured by the 
tentative shrinkage of the shock to sub-virial radii, as 
seen in the top panels, is correlated with the specific 
virial accretion rate at shock onset, which could be read 
from the red curve at the bottom panels or from the in- 
verse of the accretion time shown in Fig. [8] This trend 
is valid for the whole suite of simulated SAMBAs. 

Other quantities characterizing the tentative 
quenching and the subsequent burst are summarized 
in Figs. [6] and [71 Most of the burst characteristics vary 
gradually with the burst time % = 0.5(^2 + £3), and 
show a bimodality to early and late bursts separated at 
z ~ 1.4, associated with the shock instability. While in 
the late bursts the disc mass is Md~ 10 11 M Q , it is typi- 
cally twice as large in the early burst, where the specific 
accretion rate is higher and the preceding quenching is 
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Figure 6. SAMBA properties as a function of redshift: disc ac- 
cretion properties. Red triangles refer to the tentative quenching, 
at z\ (top panels) or z q = 0.5(^i + z%) (bottom panels). Blue 
squares refer to the burst at zj, = 0.5(z2 + 23). Vertical lines sep- 
arate early from late SAMBAs. From top to bottom: disc mass, 
fractional disc growth, accretion rate, and specific accretion rate. 



less efficient. The fraction of disc mass involved in the 
burst between ti and ta, AMd/Mj, varies from ^0.7 in 
the early bursts to ~ 0.3-0.4 in the late ones. The mean 
specific accretion rate, (AMd/At)/Md, varies from ~ 1- 
1.5 to ~ 0.4 Gyr -1 . This corresponds to a drop of the 
average accretion rate across the burst from ~ 100- 
400M Q yr _1 in the early bursts to ~ 30-40 M Q yr _1 in 
the late bursts. The average burst efficiency over the 
burst is always B h = (AM d /M d ) / (AM/M) ~ 2, while 
the peak efhciency changes from ~ 2 in the early bursts 
to 3-6 in the late bursts. 

The ratio of virial to disc mass at the end of 
the burst is bi-modal: It is high for the early bursts, 
M/M d ~ 25-50, and low for the late bursts, M/M d ~10- 
15. The high values in the early bursts is because the 
high virial accretion rate there leads to a rapid halo 
growth during the quenching phase combined with the 
failure of the subsequent burst to consume most of 
the accumulated hot gas. In the late SAMBAs, on the 
other hand, the halo growth is slow and the burst con- 
sumes all the available gas, leading to a low M/M d . The 
early bursts are thus predicted to reside in very massive 
haloes, of M ^ 10 13 M Q , while the late bursts are to be 
found in ~10 12 M Q haloes. 

The characteristics of the tentative quenching phase 
also vary gradually with t\. At the onset of all SAM- 
BAs, M - 10 12 M Q and M A ~ 10 n M Q , and we should 
scale out the apparent slight decline with time of these 
masses, which is largely an artifact of the assumed con- 



Figure 7. SAMBA properties as in Fig. \E\ disc versus virial ac- 
cretion. From top to bottom: virial mass, virial-to-disc mass ratio, 
disc-to- virial accretion ratio and efficiency B of specific-accretion 
ratio. Horizontal lines mark the fits used for SAM recipes in JS] 



stancy of metallicity with time. The typical quenching 
efficiency is B~0.3 for the early SAMBAs and B~0.4 
for the late ones. The ratio of disc to virial gas accretion 
rate, AM d /(AMf h ) varies from < 0.1 to -0.3 in the 
early and late SAMBAs respectively, separated at onset 
redshift z ~ 2.5. The tentative quenching can make the 
galaxy appear red for ~ 1 Gyr before bursting. 

The most pronounced effect of the shock instability 
occurring in the late SAMBAs is a drastic shutdown of 
disc accretion for several Gyr after the burst, starting at 
z ~ 1.2 or later. The infalling gas is halted and pushed 
back after crossing the shock during its second rapid 
expansion phase, and the fact that the halo is already 
well above M s h makes the cooling and infall time long. 
These galaxies, in haloes ^ 10 13 M Q today, will appear 
red & dead at z ^ 1, without the help of any additional 
feedback effect. On the other hand, the post-burst phase 
of the early SAMBAs is characterized by a disc specific 
accretion rate that fluctuates about the virial rate, 
1. For these central cluster galaxies today to become 
red & dead one should appeal to another quenching 
mechanism that could provide long-term maintenance. 



4 COMPARISON WITH OBSERVATIONS 
4.1 Maximum Bursting Disks 

Two relevant examples of observed "maximum 
bursters" stand out. First, rest-frame UV/optically- 
selected star forming galaxies at z ^ 2 (termed BX 
and BzK). A high-resolution adaptive-optics IR study 
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of BzK-15504 at z ~ 2.4 dGenzel et all I2006D reveals 
M d ~ 8x 10 10 M Q plus M gas ~ 4x 10 10 M Q forming stars 
at ~ 100-200 M©yr _1 . The stellar population has a typ- 
ical age of only ~0.5 Gyr, less than 20% of the Hubble 
time at that epoch. The morphology and kinematics fit 
an exponential disc of scale radius ~ 4.5 kpc rotating at 
V xo t ~ 230 km s , and a recent major merger is ruled 
out. The velocity dispersion of a ~ 60—110 kms indi- 
cates a rapid accretion onto the disc, again of ~ 0.5 Gyr. 
Several galaxies with sim ilar properties were detected i n 
a lower resolution study (|F5rster Schreiber et al.|[2006h ■ 
Other ex amples are found am ong the LIRGs at 
z ~ 0.6-0.9. iHammer et al.l (|2005[) have detected at 
least six galaxies of Mi ~ 10 11 Af Q forming stars at 
~ 120 — 200Af Q yr _1 while their morphologies and 
kinematics resemble spiral galaxies with no trace of 
major mergers. More than half of their IR-luminous 
gal axies are spir a ls. Th e compilation of SFR at z ~ 1 
by iNoeske et all f|2007h reveals some massive galaxies 
that form stars at ~ 200 Af Q yr _1 , while only ~ 7% of 
the galaxies show evidence for ongoing major mergers 
(jLotz et al.ll2006l ), and this fraction is not higher for the 
galaxies with the highest SFR (private communication). 

In summary, the observed maximum bursters show 
the following robust features: (a) the baryonic mass is 
on the order of Md ~ 10 M , (b) a substantial frac- 
tion of the baryonic mass is bursting, (c) the burst is 
short compared to the Hubble time and the cosmologi- 
cal mean accretion time, and is comparable to the virial 
crossing time, (d) the star formation occurs in an ex- 
tended configuration of a few kpc, and (e) in many cases 
there are morphological and kinematical indications for 
a thick rotating disc with no trace of major mergers. 

A distinct feature of the maximum bursts is that 
most of the stars of the M d ~ 10 Mg galaxies have 
formed over a short period of time just prior to the ob- 
served epoch, ~ 0.5 Gyr at z ^ 2 and ~ 1 Gyr at z 1. 
This is compared in Fig. [8] to several relevant charac- 
teristic timescales. The burst duration is shorter than 
the age of the universe at that time by a factor ~ 6. 
This implies that most of the stars have formed in a 
coherent starburst of the whole ~ 10 11 M© of gas in its 
late stages of collapse or assembly rather than in the 
smaller progenitors during the preceding period of hi- 
erarchical buildup. The implied suppression of SFR in 
small haloes and its possible origin will be discussed in 

m 

The burst duration is also shorter than the char- 
acteristic time for cosmological accretion onto dark 
haloes of ~ 10 12 M Q , the likely hosts of the Md ~ 
10 11 Mq maximum bursters, as given by the inverse of 
the average M/M in eq. ((!]). At z = 0.75, one has 
£acc(10 12 M Q ) ~ 7 Gyr, much longer than the indicated 
observed duration of At ~ 1 Gyr. At z — 2.4, one has 
iacc(lO 12 M0) ~ 1.6 Gyr, which is about 3 times the ob- 
served Ai~0.5Gyr, though if the relevant haloes are 
more massive, 10 13 M Q , the accretion time is ~50% 
shorter. As seen in Fig. [HI the burst durations are rather 
comparable to the shortest characteristic time of the 
system, t vn -. 

Gas-rich major mergers, which could in principle 
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Figure 8. Characteristic time scales as a function of redshift. 
The curves refer to (a) the Hubble time (blue, long dash), (b) the 
average cosmological accretion time onto M = 10 12 A/q haloes 
(red, short dash), and (c) rapid streaming from the virial radius 
to the center, R/V (green, solid). The estimated duration of the 
observed maximum starbursts are represented by black circles cir- 
cles. The predicted duration of three SAMBA bursts are marked 
by blue open squares. 



produce such rapid, massive bursts, fail to match the 
robust features (a), (d) and (e). In particular, major 
mergers are expected to destroy the discs and be as- 
sociated with perturbed clumpy morphology and kine- 
matics rather than a smooth rotating disc. This is be- 
cause the dynamical disturbance, involving tidal effects 
and violent relaxation, is expected to be pushed to an 
extreme during the ~ 100 Myr of the close passages. 
Furthermore, the induced star formation is expected to 
be concentrated in a central cusp rather than an ex- 
tended disc. It is worth noting that the mergers are 
only weakly associated with a characteristic mass, the 
nonlinear clustering scale M*, which varies rapidly with 
time — they are not associated with galaxies of any spe- 
cific fixed mass. 

The SAMBA bursts, as described by our naive 
spherical model above, have the potential of reproduc- 
ing the observed features one by one in a straightfor- 
ward way. In particular, they have the potential of keep- 
ing the disc intact while doubling its mass and puffing 
up its thickness. This is because the disc buildup, even 
if it involves filamentary and clumpy cold flows, could 
be rather uniformly stretched over the ~ 0.5-1 Gyr du- 
ration of the burst, and therefore involve only moderate 
perturbative effects on the dynamics. In addition, the 
SAMBA bursts are predicted to be specifically associ- 
ated with one characteristic baryonic mass of ~ 10 11 M Q 
at all redshifts. 

The early SAMBA bursts, which may be associ- 
ated with the observed maximum bursters at z £ 2, are 
predicted to be in M d ~2 x 1O 11 M0 discs embedded in 
rather massive haloes of M ~ 6 x 10 12 M Q . The comoving 
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number density of haloes of M > 5 x 10 12 M Q at z = 2.2 
is n ~ 10~ 4 ( /i _1 Mpc)~ 3 , consistent wit h the number 
densi ty of observed maximum bursters (|Genzel et alj 
2006)0 The predicted halo virial velocity is ~ 400 km s -1 
and the virial radius is ~ 200 kpc. The low concentration 
expec ted for such haloes at z~2.2, C~4 (Bullock et al. 
2001), implies that the circular velocity at ~ 10 kpc 
should be ~200kms — x , compatible with the observed 
bursters. The large virial radius for that epoch, with an 
assumed typical spin para meter, implies an e xponen- 
tial disc radius of ~5kpc (Bu llock et al.|[200ll ). which 
explains the observed extended discs. 

Thus, the burst properties predicted by the ideal- 
ized spherical model do surprisingly well in matching 
the observed bursts. 



4.2 Suppressed Star Formation Rate 

Despite the apparent vigor of the burst event, the most 
robust feature of the SAMBAs is the quenching. We 
note that during the SAMBA sequence of events, from 
the onset of the shock till the end of the burst about a 
Hubble time later, each galaxy spends more time in a 
phase of suppressed accretion rate rather than in a rapid 
accretion mode. Figure [H] shows a crude estimate of the 
predicted distribution of disc accretion rate for galaxies 
in the stellar mass bin (0.5-2) x 10 11 M Q near z~ 2. This 
crude estimate is based on the fraction of time spent at 
each accretion rate in the simulated SAMBA shown in 
Fig. [31 assuming that it is a representative SAMBA at 
that epoch. While a fraction of the galaxies are expected 
to show a high accretion rate into the disc at about twice 
the average virial level, more than half the galaxies are 
predicted to present a strongly suppressed disc accretion 
rate, an order of magnitude below the virial accretion 
rate or less. The associated distribution of SFR may be 
qualitatively similar. 

Preliminary indications for a substantial population 
of ~ 10 11 Mq galaxies with strongly suppressed SFR 
already exist. A bi modality in the SFR distribution is 
apparent at z 1 in lNoeske et al.l ()2007l Fig. 1), where 
the typical SFR is > lOAf yr _1 while about half the 
galaxies h ave an upper limi t of SFR< lMoyr -1 . At 
z ~ 2-2.7, iKriek et alJ (l200ll report that a significant 
fraction of the ~ 10 n MQ galaxies have surprisingly low 
upper limits for the specific SFR at < 0.1 Gyr -1 , while 
the typical is ~1 Gyr" 1 . 

The SAMBAs predict that the SFR histories of 
galaxies should show two peaks of SFR separated by 
a couple of Gyr, but this might be hard to detect ob- 
servationally. 



2 This estimate is for the WMAP3 parameters of the ACDM 
cos mology, H m = 0.2 4, Q\ = 0.76, h = 0.73, and <rg = 
0.74 jSpergel et al . 2006). A similar number density is obtained for 
haloes twice as massive if one adopts instead f2 m = 0.3, Q\ = 0.7, 
h = 0.7, and erg = 0.9. 
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Figure 9. Predicted distribution of log disc accretion rate for 
galaxies in the disc-mass bin (0.5-2) X 10 11 Mq near z ~ 2. The 
vertical green line marks the mean virial accretion rate. About 
half the galaxies are expected to show a high accretion rate into 
the disc, while the other half is predicted to suffer a strongly 
suppressed accretion rate. 

4.3 Downsizing 

The SAMBA phenomenon provides a simple explana- 
tion for the observed "downsizing" of elliptical galax- 
ies. An "archaeological" analysis of stellar ages reveals 
that today's more massive ellipticals ha ve formed their 
stars earlier and over a shorter time span ([Thomas et alj 
120051) . From their SFR histories as a function of look- 
back time one learns that elliptical galaxies of present- 
day stellar mass ~ 10 11 — 10 12 M Q have formed most of 
their stars at z~2-3 in a maximum burst over 1-0.5 Gyr 
with Md/Md ~ 2 (or z ~ 1.2-1.6 in lower density envi- 
ronments). The progenitors of today's ellipticals, prior 
to their subsequent growth via dry mergers, seem to 
be consistent with the maximum bursters seen at high 
redshift, and with the SAMBA predictions. 

The constancy of M s h with redshift, which implies 
that all galaxies burst when they are 1O 11 M0 and 
quench immediately after automatically implies down- 
sizing in both the peak of SFR and the subsequent 
quenching. Galaxies of larger stellar mass today are 
likely to be embedded in more massive haloes, which 
in turn have crossed the threshold mass for the onset 
of SAMBA at an earlier time. They bursted earlier, for 
a shorter duration, and shut down earlier accordingly. 
This is demo nstrated in a cosmolo gical semi-analytic 
simulation bv lCattaneo et aD (|2007f ). 

The SAMBAs predict yet another kind of down- 
sizing phenomenon, concerning halo mass at burst: the 
early bursts at z > 1.4 are predicted to reside in mas- 
sive haloes only slightly smaller than Al ~ 10 13 M Q , 
while the late bursts are to be found in smaller haloes, 
M~10 12 M ra . 



5 RECIPES FOR SEMI- ANALYTIC MODELS 

In order to study the potential observable implica- 
tions of the SAMBA phenomenon in more detail, one 
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Figure 10. Disc growth rate in the SAMBA recipe for SAM 
(blue, solid) versus the smooth virial baryonic accretion rate (red, 
long dash) and the SAM standard (gree n, short dash) b a sed o n 
the evolution of the "cooling radius" a la I White fc Frenk] l|l99ll ). 
Top: an early SAMBA {z\ ~ 4.8, M ~ 10 14 M Q ). In this case, 
long-term quenching has to originate from another process. Bot- 
tom: a late SAMBA (z 1 ~ 1.8, M ~ 4 X 10 12 M Q ), were long- 
term quenching is a natural outcome of uniform accretion. 



can translate the predictions from the spherical model 
into a recipe for gas cooling and accretion rate onto 
the disc. This recipe is to be incorporated in Semi- 
Analytic Models (SAM), which follow the baryonic 
physics in given dark-halo merger trees. The SAMBA 
recipe should replace the st andard recipe adopted in 
SAMs rtWhite fe Frenldll99l . which assumes that the 
halo gas interior to a "cooling radius" is accreted onto 
the disc. The cooling radius within a (spherical) halo 
of a given mass at a given time is traditionally deter- 
mined by the position where the cooling time equals 
a certain dynamical time, e.g., the Hubble time. The 
SAMBA phenomenon suggests a deterministic modifi- 
cation of the disc accretion rate, based on our spherical 
modeling, with no free parameters. However, an optimal 
match to observations may suggest a certain fine-tuning 
of the model parameters, compensating for deviations 
from spherical symmetry and uniform accretion. 



5.1 Following Halo Evolution in SAM 

The proposed SAMBA recipe for the accretion rate onto 
the disc is specified in £)B1I In summary, given the mass 
growth of each halo in the merger tree, M(t), and its 
derivative M, the SAMBA onset time t\ (redshift z\) is 
identified by M(t x ) = M x = 1O 12 M , and the SAMBA 
is classified as "early type" or "late type" according 
to whether z\ > 2.5 or z\ ^ 2.5 respectively. Then 
the times ti and 13, the beginning and the end of the 
burst, are computed using the fits shown in Fig.O The 
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Figure 11. Disc growth in the SAMBA recipe for FHM (blue, 
solid) versus the virial baryonic accretion rate (red, long dash) 
and the standard based on "cooling radius" (green, short dash). 
Top: at z = 2. Bottom: at z = 1. Quenching in haloes ~ 10 13 Mq 
and larger has to be due to another process. 



algorithm then specifies the growth rates of cold and 
hot gas mass during each of the four SAMBA phases 
defined by the above times. The numerical factors used 
for the tentative quenching and for the burst are read 
from straightforward fits to AM d /(AM/ b ) in Fig. [3 

The disc mass growth according to the proposed 
SAMBA recipe for SAM is shown in Fig. [TO] for typical 
early and late SAMBAs. It is compared to the standard 
disc growth rate following the evolution of a "cooling 
radius" . The latter has a characteristic drop when the 
cooling radius first becomes smaller than the virial ra- 
dius. The cooling-radius recipe can be somewhat im- 
proved once the Hubble time used in its definition is 
replaced by i V i r = 0.18iHubbie- Yet, it fails to capture 
the depth of the tentative quenching, the burst, and the 
long-term quenching of late SAMBAs. 

Recall that the early SAMBAs, ending up in M ^ 
10 13 M Q haloes today, should be quenched after t 3 by 
another mechanism, such as clumpy accretion or AGN 
feedback. If such a process is not explicitly included in 
the SAM, one can simply set Md = in the po st-burst 
phase for all haloes (as in lCattaneo et ai1l2006[ ). 



5.2 SFR as a function of Halo Mass: FHM 

An alternative, quick and crude way to learn about the 
implications of SAMBAs is by providing a recipe for the 
SFR as a function of halo mass and redshift, to be in- 
corporated in a given halo population in a partial SAM 
while bypassing the detailed baryonic processes of cool- 
ing, star formati on, and feedback (For mation History 
Modeling, FHM. ICroton fe Dekell [20071 ). The fact that 
the SAMBA phenomenon is driven by halo mass allows 
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us to provide a concrete recipe of this sort, specified in 

In summary, one starts with M\ — 10 12 M Q at ti, 
and computes the other SAMBA characteristic halo 
masses Mi and M3, corresponding to the characteris- 
tic times ti and £3, using fits to our simulation results. 
Given the virial accretion rate in the provided merger 
tree, M, the disc accretion rate, which could be inter- 
preted as an upper limit for the SFR, is given in each 
of the mass ranges, corresponding to the four SAMBA 
phases. The numerical values are read from crude fits to 
the simulation results for AMd/(AM/b), as shown in 
Fig. [71 When implementing the recipe, one can either 
use the actual mass growth rate of the haloes drawn 
from merger-tree realizations, or the average virial ac- 
cretion rate as provided in eq. ©. 

The SAMBA recipe for disc accretion rate is shown 
in Fig. [TTl compared to the standard rate based on a 
"cooling radius". Again, the cooling-radius recipe can 
be somewhat improved once the Hubble time is replaced 
by t v ir, but it would still fail to recover the burst and 
the deep quenching features of the SAMBAs. 

The main feature of the SAMBA FHM is the 
quenching of SFR above a threshold halo mass, M s h ~ 
10 12 M(T), as already implemented in SAM simulations 
(jCattaneo et all 120061 : ICroton et aI1l2006h . The inter- 
esting new feature predicted by the model is the high 
SFR in a second range of masses above M 8 h. At low 
redshifts, this is a slight increase of the effective thresh- 
old mass for quenching to ~ 2 x 10 12 M Q . At z ^ 1.3, 
the high SFR is predicted to extend to haloes as large 
as ~ 7 x 10 12 M Q . This should allow the appearance of 
massive blue galaxies at high redshifts. In the special 
range 1.24 < z < 1.53, the high SFR is expected over 
an especially broad range of masses, 1.5-7 x 10 12 Af Q 
(not shown in Fig. ITT)) . Given that 10 12 M© is the typ- 
ical mass for haloes that form at this redshift range, 
this should lead to a peak in the history of cosmological 
SFR density (the "Madau" plot). 



6 DISCUSSION AND CONCLUSION 

A straightforward spherical model of uniform gas ac- 
cretion onto haloes predicts a robust, rich sequence of 
events in the accretion onto the central discs, SAMBA, 
which could be the major driver of star-formation his- 
tory in galaxies. An initial cold accretion phase ends 
abruptly upon the birth of a shock in the inner regions 
of a halo once its mass exceeds a threshold mass of 
~ 10 12 M© and the disc contains ~ 1O U M . The infall 
of gas through the expanding shock creates a period of 
tentative quenching, where the accretion onto the disc 
is suppressed for a couple of Gyr. The accumulating hot 
gas joins new incoming gas in a big crunch, leading to a 
massive, rapid, burst involving ~ lO n M0 of gas in a few 
hundred Myr. The subsequent expansion of the shock to 
the virial radius results in a long-term shutdown, which 
is especially efficient if the burst occurs after z ~ 1.4, 
corresponding to today's haloes of ^ 1O 13 M0. The pre- 
dicted SAMBA bursts seem to match the observed max- 
imum bursters of Md ~ 10 11 M Q at z £ 1 and at z 2, 



as well as the indicated population of Md ~ 10 11 Mq 
galaxies with suppressed SFR. This rich phenomenon is 
a generic consequence of uniform accretion, before in- 
cluding the effects of major mergers or AGN feedback. 

The implications of the SAMBA phenomenon 
should be taken with a grain of salt because the rel- 
evance of the spherical uniform accretion model to 
the actual SFR history in real galaxies is yet to be 
demonstrated. This is an immediate challenge for cos- 
mological simulations. The actual accretion could be 
rather filamentary and clumpy to begin with, and one 
should verify the extent to which the SAMBA behav- 
ior remains valid under such conditions. One could 
be encouraged though by the success of the spheri- 
cal model in pr edicting the threshold mass for virial 
shock heating dKeres et all 120051 : iDekel fc Birnboim! 
I2006t iBirnboim et all I2007D . First clues for natural 
long-term quenching of the sort predicted here in a 
~ 10 12 M Q halo may be seen in high- resolution simula- 
tions of galaxies with qu iet histories (|Naab et al.l 120071 : 
iLibeskind & Dekell l2007t) . The following discussion of 
possible implications is pursued under the tentative as- 
sumption that other robust SAMBA features will also 
be reproduced in cosmological simulations of the forma- 
tion of such galaxies. The apparent match of the spheri- 
cal predictions with several different observations makes 
this discussion worthwhile. 

Using the fact that the accretion onto the disc is a 
straightforward function of halo mass, we propose a sim- 
ple accretion recipe for SAMs, to replace the standard 
recipe based on a "cooling radius" . A realistic study of 
the SFR that is associated with this galaxy buildup, 
both in the case of smooth accretion of minor mergers 
and under major mergers and feedback effects, is neces- 
sary for a reliable comparison with observations. Still, 
the naive interpretation of the disc accretion rate as the 
maximum possible SFR provides a preliminary indica- 
tion for a surprisingly good match to observations. 

An obvious necessary condition for SAMBAs 
(which is actually automatically implied by the ob- 
served maximum bursts independent of any model) is 
that a substantial fraction of the accreting material onto 
haloes of ~ 10 12 M Q is gaseous rather than stellar. Star 
formation could have been suppressed in haloes smaller 
than a certain threshold due to different processes, such 
as stellar feedback or a density threshold for star for- 
mation. Given such a threshold mass, we can estimate 
the gas fraction in the accretion by evaluating the mass 
in all accreting haloes bel ow that threshold comp ared 
to the total accreting mass dNeistein fc De kcl 2007). We 
find using the EPS formalism specified in <TA]that about 
one half of the accreting mass onto a halo of 10 12 M Q is 
in haloes smaller than ~ 10 11 AI Q (or "smooth" accre- 
tion), independent of time. This should therefore be the 
minimum halo mass for star formation in order to have 
a gas-dominated accretion. If the suppression is due to 
supernova feedback, the minimum mass for star forma- 
tion corresponds to haloes of virial velocity ~ 120 kms -1 
(IDekel fc Silldll986blDekel fc Birnboimll2006l Fig. 2), for 
which the accretion is almost all gaseous at z = 0, it is 
about 90% at z ~ 0.7, and about 67% at z ~ 2.4. This 
hints that supernova feedback could be the suppression 
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mechanism allowing the accretion onto ~ M s h haloes to 
be predominantly gaseous at all relevant redshifts. 

The dramatic variations in the disc accretion rate 
imply that each galaxy is likely to evolve through sev- 
eral subsequent phases of very high and very low SFR. 
This means that the expected overall trend of evolution 
from the blue cloud to the red sequence of galaxies is 
not a one-way track. In particular, some of the red & 
dead galaxies observed at high redshift may actually be- 
come blue again at later redshifts. The lesson is that the 
interpretation of the observed evolution of the blue and 
red luminosity functions may involve non-trivial book- 
keeping. 

The late SAMBAs offer a very efficient post-burst 
shutdown that lasts for 6-7 Gyr. This can naturally give 
rise to red & dead galaxies in today's small-group haloes 
of < 10 13 M(D, perhaps SO or elliptical galaxies (as si m- 
ulated bv lNaab et atll2007t ILibeskind fc Dekelll2007ft . 

In the early SAMBAs, the accretion rate into the 
disc is suppressed compared to the virial rate even dur- 
ing the burst phase, but the post-burst accretion rate is 
not negligible. For these central galaxies in today's clus- 
ter haloes of > 10 13 M Q to end up as large red & dead 
ellipticals one needs a quenching mechanism beyond the 
one provided by spherica l accretion. We show elsewhere 
(|Birnboim fe Dekell2007l ) that this could be provided by 
clumpy accretion, in which gas clumps of 10 7 — 1O 9 M 
transfer their gravitational energy via drag and local 
shocks into heating and puffing up of the inner-halo gas 
of ^ 10 13 M Q haloes. 

Alternatively, a fashionable scenario suggests that 
the required quenching is provided by A GN feedback 
(jCroton et al.N2006t ICattaneo et al.ll2006l) . It has been 
noticed that the shock heating of the gas in haloes 
> 10 12 M Q enables the r equired coupling of the AGN 
energy with the halo gas (| Dekel fc Birnb oim 2006|, but 
we now realize that the SAMBA burst may be the ac- 
tual trigger for AGN activity just above this critical 
mass. The SAMBAs may thus be responsible for both, 
the energy source and the coupling mechanism. 

If the SAMBAs are responsible for triggering AGN 
activity, recall that they do so only in galaxies of 
~ 10 11 M Q , a couple of Gyr after their haloes became 
larger than ~ 10 12 M Q . This may explain the other- 
wise surprising result that the bright quasars (near L* 
of their Schechter luminosity function) s eem to reside 
at all redshifts in haloes of ~2 x 10 12 M Q (|Croom et al.l 
12005ft . A straightforward calculation using the Press- 
Schcchtcr formalism shows that the evolution of num- 
ber density of quasars at z > 1 can be reproduced by 
assuming quasar onset in every halo when it first be- 
comes ~2x 10 12 Af Q and allowing a fixed short lifetime 
for each quasar. 

The SAMBA events offer a complementary alterna- 
tive to wet major mergers in producing rapid, massive 
bursts, with the important difference that it can leave 
the disc intact. The associated shutdown leads to red 
& dead galaxies in haloes above the threshold mass. 
However, while major mergers also automatically pro- 
duce spheroids, what could be the mechanism respon- 
sible for spheroid formation in association with the old, 
gas-poor stellar population according the halo shutdown 



scenario? First, the SAMBA burst provides the gas nec- 
essary for bar instability in the disc, which then trans- 
fers angular momentum to the outer disc and halo and 
produces a bulge. Second, the SAMBA quenching makes 
any subsequent mergers dryer, so they could lead to the 
large, boxy, non-rotating elliptic als, which cannot be re - 
produced by wet major mergers (jNovak fc et al.][2007l ). 
Third, the minor mergers associated with the contin- 
uous smooth accretion of SAMBA could themselves 
lead to the formation of a spheroid ( Naab et all 20071 : 
iBourneaud fc et al. J 120071: ILibeskind fc Dekelll2007[) . 

We argued in iDekel fc Birnboirq (|2006l ) that discs 
could be built and form stars by cold streams, which 
could persist as narrow filaments even within the shock 
heated medium slightly above M s h, especially at high 
redshifts. Streaming at the virial velocity, these flows 
could indeed give rise to rapid starbursts, on a time 
scale comparable to the virial crossing time, but it is not 
clear that these flows could bring in mass in a rate that 
is much more efficient than the cosmological accretion 
rate. The SAMBAs provide such more efficient bursts 
in a natural way. 

All the above is yet to be confirmed in cosmological 
simulations. It would be interesting to learn how the 
SAMBA phenomenon manifests itself when the spheri- 
cal accretion of minor mergers is accompanied by fila- 
mentary streams and major mergers. 
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the cosmological growth factor of linear density fluc- 
tuations normalized to D(to) = 1 (see below), and 
S c ~ 1.68 (ignoring small variations due to the actual 
cosmology used) . Mass enters via the monotonically de- 
creasing function a(M), the standard deviation of the 
initial density fluctuations, as derived from their power 
spectrum, smoothed on a scale that corresponds to mass 
M and linearly extrapolated to to (see below). Accord- 
ing to EPS, a halo of mass Mq at luq has the follow- 
ing average number of progenitors of mass in the range 
(M, M + dM) at the earlier time uj + Auj: 



diV=- 



1 



Alu 



/2¥ (Aa 2 ) 3 / 2 



exp 



2Aa 2 



da 2 



dM 



^dM.(Al) 



where Alu = uj — uj and Act 2 = a 2 (M) — a 2 (M ). 



Al Accretion Rate onto M at t: a Small Time 
Step 

We wish to compute the main progenitor history, i.e., 
the conditional average M(lu\luq, Mq). The main pro- 
genitor is guaranteed to also be the most massive pro- 
genitor when Aw is small enough, and this allows us to 
compute M by evaluating M(u>o + Alu\luq, Mo) for small 
Auj. 

For M ^ M /2, the probability density 
P(M, lu\Mq, luq) is the same as the total progenitor 
distribution dN/dM given by equation IA1I For M < 
Mo/2, however, the only condition is P ^ dN/dM, and 
it is not sufficient to constrain P. As a first approxima- 
tion we assume that the main progenitor always has a 
mass M ^ Mo/2. Then, the average mass of the main 
progenitor is 



M(w + Aoj\ui ,M q 



MdM 



= M 



1 - erf 




(A2) 



where a q — a(Mo/q) with q = 2, and <tq — a(Mo). The 
rate of change is then 



dM 
dto 



lim 



M (ujq + Aw|w , M ) - M 



—Mq lim — — erf 

Aw^o Alu 




(A3) 



Using the limit erf(x) 



1 dM 

M duj 



2/tt 



► 2x/y / 7r when x 

1/2 



we get 

(A4) 



APPENDIX A: MAIN PROGENITOR HISTORY 



Following iNeistein et al l (120061) . we use the EPS 



formalism (|Lacev fc Coldll993l ) to derive an analytic es- 
timate for the main progenitor history. We track time 
backwards via the variable u> = 8 C /D(t), where D{t) is 



Thus, the mean specific accretion rate onto haloes of 
mass M at t is 



M , 

A? 

with 



-s(M)ib(t) 



(A5) 
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s{M) 



w(t) 



2/tt 

a 2 (M/q)-a 2 (M) 

ScD 
D D ' 



1/2 



(A6) 



(A7) 



When Aoj is not negligibly small, the probability for 
M < Mp/2 is not neg l igible , as it grows in proportion 
to Auj. iNeistein et aTT (|2006l ) showed that this effect is 
rather small and can be bounded by well defined limits 
on q, though the exact value of q is not fully specified by 
the EPS formalism. They showed that eq. (|A5[) is still 
a valid approximation, but with q slightly larger than 
2, and limited to the range 2-2.3 for a flat ACDM cos- 
mology with 0.1 < il m < 0.9 and 10 8 < M < 1O 15 M . 
We adopt q = 2.2 for all practical purposes. 

As a sanity check, note that the obtained specific 
accretion rate at a fixed mass is self-similar in time, 
consistent with the fact that the only relevant mass 
scale at any given time is the non-linear clustering 
mass of the Press-Schechter formalism, M*, defined by 
cr[M*(i)] = 8 C /D(t). This is because the time depen- 
dence in eq. (|A5|) is via Co = —cr{M*(t)]D/ D. Once 
M/M is divided by the universal growth rate D/D, the 
time dependence enters only through M*(t). Any viable 
approximation for the cosmological accretion rate must 
obey a self-similarity of this sort. 

We also note in eq. ^ that at a fixed t, M/M — » oo 
when M — > oo, but recall that haloes of M ^ M* are 
rare, such that most of the accreted mass is not onto 
arbitrarily massive haloes. 

A2 Mass Growth History of M at t 

The mass growth history of the main progenitor of a 
halo of mass Mq at to is obtained by integrating eq. (j A4[) 
over du> and over dM in the two sides of the equation: 



u-uq = F(M)-F(M ), F{M) 



dm 



(A8) 



l M ms(m) 

The desired average mass of the main progenitor at t 
can be extracted from the above equation, 

M(t\M ) = F->-u, + F(M )}, w = (A9) 

where F^ 1 is the inverse function of F, to be evaluated 
either explicitly or via numerical interpolation. 

A3 A Power-Law Power Spectrum 

A simple, fully self-similar example is provided by a 
power-law power spectrum, oc k n , for which 



<r = cr R M~ a , a=(n + 3)/6. 



(A10) 



Here or is a normalization constant, and M is measured 
in units of Mr, the average mass originally contained 
in a top-hat sphere of comoving radius R (e.g. as for 
R = 8 ft. Mpc). In this case, the average accretion rate 
is given by eq. (|A5|) with 



s(M) = rM a , t= — 



2/tt 

2a 



1/2 



(All) 



Then, 



F(M) = — (M~ a - Mq 
ar 



(A12) 



and the explicit solution for the average main- 
progenitor mass is 



M(t\M ) = [M - a + a T {uo - u T 1/a 



(A13) 



(with to = S c /D(t), M in units of Mr, a nd q ~ 2.2). 

Eq. (|A5|) with s(M) from eq. (|A1 1|) approximates 
the accretion rate for a general power spectrum once a 
is the local slope a(M) = —d\na/dlnM. For n = —2.1, 
appropriate for M~10 12 M Q in the ACDM cosmology, 
we have a ~ 0.15. This is indeed similar to the power 
index of the mass dep endence from N-body simulations 
(|Wechsler et alJl2002L Fig. 12). 

The time dependence can be approximated in stan- 
dard ACDM by eq. l[A29]l . Using a crude power-law ap- 
proximation for the growth of the Press-Schechter mass 
out to z^2, M*(z)/M* — (1 + z)~ 5 , we end up with 
a practical approximation for the accretion rate into 
haloes in the vicinity of M~ 10 12 M Q , 



M 
~M 



0.04 Gyr~ 1 {l + z) 2 - 25 M- 



12 : 



(A14) 



which is accurate to better than 10% over the range 
^ z ^ 2. 



A4 ACDM Power Spectrum 

For the standard ACDM power spectrum we 
adopt P(k) = kT 2 (k) with the transfer function 
(|Bardeen et al.lll986l ) 



T(k)- 



ln(l + 2.34 k) 



2.34 k 

x[1 + 3.89k+(16.1k) 2 +(5.46k) 3 + (6.71k) 



(A15) 

41-1/4 



Here k = k/T, with the wave number k in units 
of (fe~ 1 Mpc) , and with the shape parameter T 
(|Sugivamalll995h 



r = ^ m /iexp -n h (i + V2h/n„ 



(A16) 



This power spectrum can be used to derive cr(M) and 
the corresponding s(M) with q = 2.2 in eq. (J5J). This can 
then be substituted in eq. (|A5[) for an explicit estimate 
of the average accretion rate onto haloes of mass M at 
t. It can then be used in eq. I|A8[) for computing F(M), 
evaluating its inverse function numerically, and using 
it in eq. ((4]) to obtain the average history of the main 
progenitor, M(t\M ). 

Alternatively, following Ivan den Boschl (|2002T> . the 
rms fluctuation corresponding the ACDM power spec- 
trum is approximated (to better 1% in the range 10 6 < 
M < 10 16 Af o ) by 

^ = ^«1, (AIT) 
M c = 9.73 x 10 n M Q ho$ (fim/0.3) (r/0.2)" 3 , (A18) 
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and with the fitting function 
*(x) = 64.09 (1 + 1.074 a; - 3 - 1.581 x 0A 
+0.954x - 5 - 0.185 a; - 6 )" 10 . 



(A19) 



Substituting this analytic approximation in 
eq. (|A8p . numerically in v erting F(M), and using it 
m cq. (@]), iNeistein etail (|2006h derived the following 
approximation for the average main progenitor history: 



M(t\M ) 



(W-Wo) J rF q [ — 



, (A20) 



with q — 2.2 and the fitting function 

F2. 2 {u) = -6.92 x 10~ 5 (lnu) 4 + 5.0 x 10~ 3 (lnu) 3 (A21) 
+8.64x 10~ 2 (lnu) 2 - 12.66 Qnu) . 

This fitting function is accurate to better than 1% over 
the range 10 6 < M < 1O 15 M for the standard ACDM 
cosmology. 



A5 In Cosmological N-body Simulations 

IWechsler et al.l (|2002T ) analyzed an N-body simulation 
of ACDM for the average mass-growth history of the 
main progenitor of a halo of mass Mo at to . They found 
a good fit with 

M(t\M ) = M e~ 2a * z , 



a c = 0.28 



(A22) 
(A23) 



This function is a good practical approximation to the 
EPS estimate specified above, as long as t is the present 
time, z = 0. 



A6 Useful Relations for Linear Fluctuations 

Followin EtMo fc White! (|2002f) . in a flat universe 

9(z) 1 



D(z) 



5 (0) (1 + z) 



(A24) 



(A25) 



with a negligible error on the order of Q\/70 in g(z). 
Thus 



H(z) (l + ^ g ( z )n-W( z )[n m (z)-lfj . (A26) 



D ^ 
D ~ 

Here 

H(z)=H E(z), n m (z) = nm f 2( \ Z)3 , (A27) 

E 2 (z) = l-n m0 + n rn0 (l + z) 3 . (A28) 

For ACDM with fi m = 0.3 and f2 A = 0.7, a practical 
approximation is 



D 
D 



0.0385 (l + z) 3/2 ho.TGyr" 1 , 



(A29) 



which is accurate to better than 1% at z 1 and better 
than 10% at z > 0. 

The fluctuation variance is derived in the general 
case by integrating the power spectrum, 



1 

2n Jo 



a 2 (M) = —l dkk 2 P{k)W 2 (kR) , 



(A30) 



with M = (ATr/3)poR i , and with the Fourier transform 
of the real-space top-hat window function or radius R 

W{x) =3{sinx - xcosx)/x 3 . (A31) 



APPENDIX B: RECIPES FOR SAMS 

Bl Following Halo Evolution in SAM 

Given the mass growth of each halo in the merger tree, 
M(t), the SAMBA onset time h (redshift zx) is identi- 
fied by M(ti) = Mi = 10 12 M Q . The SAMBA is classi- 
fied as "early type" or "late type" according to whether 
Z\ > 2.5 or z\ ^ 2.5 respectively. Then the beginning 
and the end of the burst, as read from Fig. [51 are 



ti = 



{! 



1.80t x , 
36t x , 



zi > 2.5 
zi ^ 2.5 



h = 1-18*2 



(Bl) 



This defines the four phases of gas cooling. 

The growth rates of cold and hot gas mass during 
each phase can be approximated as follows: 

1. Cold- flow phase, t ^ t\, 



M d = f h M , Mhot = . 

2. Tentative-quenching phase, t\ < t ^ t 2 , 

■ _ AM 12 r 0.08, z x > 2.5 
Md_/b A*77 X l0.3, 2l <2.5 ' 

M hot = fhM - M d , 



(B2) 



(B3) 
(B4) 



where AM {j = M(tj)-M(ti) and = The nu- 

merical factors for early and late SAMBAs, both for the 
tentative quenching and for the burst below, are read 
from straightforward fits to AMd/(AAf/ D ) in Fig. [71 

3. Burst phase, t 2 < * ^ *3, 



, v , AM 23 
M d = / b -r- — x 

^*23 
M hot2 



/ 0.5, z x > 2.5 
1 1.6, zi s$ 2.5 



M, 



hot 



A* 



2:1 



(B5) 



(B6) 



where Mhot 2 is t ne accumulated hot gas mass by t 2 , 
Mhot2 = AAM12 — Mdisc,2A*i2. The cold accretion is 
more efficient in the late bursts because of the tentative 
shock instability due to the lower virial accretion rate 
at late times. 



4. Post-burst phase, t > £3, 

0.12, zx > 2.5 



M d = hM x { q" 
Mhot = fhM - M d 



zx < 2.5 



(B7) 
(B8) 
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The numerical factor for early SAMBAs is read from the 
fit M/Md ~ 4//b at to the results shown in Fig. [3 com- 
bined with the fact, seen in Fig.rj] that during the post- 
burst phase M d /M d ~ M/M and M/M d ~ const. Re- 
call that the early SAMBAs, ending up in M > 1O 13 M 
haloes today, should be quenched after £3 by another 
mechanism, such as clumpy accretion or AGN feedback. 
If such a process is not explicitly included in the SAM, 
one can simply set Md = in the po st-burst phase for 
all haloes (as in lCattaneo et al.ll2006f) . 

The disc mass growth according to the proposed 
SAMBA recipe for SAM is shown in Fig.[Tni in compar- 
ison with the standard disc growth following the evolu- 
tion of a "cooling radius" . 

B2 SFR as a function of Halo Mass: FHM 

An alternative, quick and crude way to learn about the 
implications of SAMBAs is by providing a recipe for 
the SFR as a function of halo mass and redshift, to 
be incorporated in a given halo population in a SAM 
while bypassing the detailed baryonic processes of cool- 
ing, star formati on, and feedback (For mation History 
Modeling, FHM. ICroton fc Dekell [20071) . The fact that 
the SAMBA phenomenon is driven by halo mass allows 
us to provide a concrete recipe of this sort. 

The SAMBA characteristic halo masses Mi, M2 
and M3, corresponding to the characteristic times t\, 
ti and £3, can be approximated based on our simula- 



tions by 

Ah = 1O 12 M , (B9) 



The constant value o f Mi is as predicted by 
iDekel fc Birnboiml (|2006f) for a realistic gradual increase 
of metallicity with time. Then M2 and M3 are obtained 
from the fits shown in Fig. IBll Notice the interesting 
prediction that in the range 1.24 < z ^ 1.53 M3 is much 
larger than Mi , giving rise to a broad range of massive 
haloes with high SFR. 

Given the virial accretion rate in the provided 
merger tree, M, the disc accretion rate, which could 
be interpreted as an upper limit for the SFR, is approx- 
imated in the different mass ranges by: 

1. Cold-flow range, M < Mi, 

M d = / b M . (B12) 

2. Tentative-quenching range, Mi < M ^ M2, 

Md = / b A/x{™ 8 < l>{* . (B13) 

3. Burst range, M 2 < M ^ M 3 , 

M d = /bMx{°J Z z l\ A 4 . (B14) 

The numerical values are read from crude fits to the sim- 
ulation results for AMd/(AM/b), as shown in Fig. [7] 
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Figure Bl. The characteristic SAMBA masses relative to Mi ~ 
10 12 Mq. Open red triangles refer to z\. Solid blue squares refer 
to M2/M1 at Z2- Solid green circles refer to M3/M1 at 23. The 
lines mark the fits used for the FHM recipe. 



4. Post-burst range, M > M 3 , 

»> t »> f 0.12, M > MJz) /TJ1K , 
M d = / b Mx| 0) Af ^ Mq(z) ' (B15) 

M q (z) ~ 1O 13 M exp(-0.75 z) . (B16) 

The upper-limit for effective quenching, M q (z), is the 
mass of the intermediate SAMBA separating late from 
early SAMBAs, which eventually ends up as Al = 
10 13 M Q today. The above estimate is based on the 
fit to simulations, eq. ()A22|) . Note that once z > 1.4, 
M3 > M q (z), so in this redshift range M > M3 auto- 
matically implies M > M q . 

The numerical values in the quenching and burst 
phases are read from the crude fits to the simulation 
results shown in Fig. [71 in most of which the virial 
mass grows according to the average cosmological rate. 
When implementing the above recipe, one can either use 
the actual mass growth rate of the haloes drawn from 
merger-tree realizations, or the average virial accretion 
rate as provided in eq. ([3]). 

The main feature of the model is the quench- 
ing of SFR above a threshold halo mass, M s h ~ 
10 12 M(7 1 , as already implemented in SAM simulations 
(|Cattaneo et all l2006t ICroton et aLll2006fh The inter- 
esting new feature predicted by the model is the high 
SFR in a second range of masses above M s h- At low red- 
shifts, this is a slight effective increase of the threshold 
mass to -2 x 10 12 M Q . At z ^ 1.3, the high SFR is pre- 
dicted to extend to haloes as large as ~ 7x 10 12 M Q . This 
should allow the appearance of massive blue galaxies at 
high rcdshifts. In the special range 1.24 < z < 1.53, the 
high SFR is expected over an especially broad range 
of masses. Given that 10 12 M Q is the typical mass for 
haloes that form at this redshift range, this should lead 
to a peak in the history of cosmological SFR density 
(the "Madau" plot). 
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B3 Fine-tuning of the SAMBA Recipes 

While the above recipes are our best guess based on 
the spherical model, with no free parameters, a better 
match to observations may require modest fine-tuning 
of model parameters, e.g., allowing for certain uncer- 
tainties in the assumptions made in the computations, 
as well as compensating for deviations from spherical 
symmetry and uniform accretion. 

The uncertainty in the value of the critical mass 
for shock heating, M s h, can easily accommodate a fac- 
to r of two and some redshift dependence. As computed 
in lDekel fc Birnboim! (|2006[ ), it is a function of certain 
physical parameters, in particular redshift z and the 
average metallicity Z and gas fraction /b in the inner 
halo, which also vary with z. The estimates of M s h as 
a function of z are shown in their figure 2. A practical 
approximation based on their equation 34 is 

,,2,, zSffn'A /lna + a,;^' 4 



"* aw, ^^m(wm) ■ <Bi7) 

where Z . 2 = Z/( 2Z m ) and f hn n* = fb/0 -05. In this 
formula, following iDekel fc Birnboiml (|2006f) . the vari- 
ation of mean metallicity with redshift can be crudely 
approximated by 

log^- = -0.17 z. (B18) 
Zo 

Following iBullock et al.l (|200ll ) , the average halo con- 
centration near M~10 12 Mq is evolving as 

C = C (l + z)- 1 , C ~13. (B19) 

With this choice, the critical mass is M s h ~ 10 12 M Q 
and is only a weak function of redshift. It becomes even 
weaker if /b gets slightly higher at high z. One can see 
that a factor of 4 change in the unknown Zq results in 
a factor of ~ 2 in M s h . 

Another uncertainty, which is explicit in the FHM 
recipe, is in the actual SFR that is associated with the 
maximum possible rate provided by the disc accretion 
rate. This can be modeled as a free ratio, smaller than 
unity, between the two. 



